A New Quantum Phase Transition in the Coupled Quantum Dots System 
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We study two quantum dots in the limit of strong dot-lead coupling and weak dot-dot tunneling. 
The model maps on Ising-coupled Kondo impurities. We argue that a new quantum critical fixed 
point exists at an intermediate value of the mutual capacitance, supporting non-Fermi liquid be- 
haviour. We construct the total conductance across the double dot structure. It exhibits a strongly 
peaked behaviour as a function of the mutual capacitance, gate voltage, and temperature. 
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Electron tunneling through quantum dots is funda- 
mentally affected by intriguing many body effects. The 
Coulomb interaction imposes a prohibitive energy cost 
Ec on the transfer of electrons, known as Coulomb block- 
ade |l],H ■ Fine tuning of the gate voltage Vq is required 
to reinstate charge flow, manifesting itself in sharp con- 
ductance peaks as a function of Vq. 

Remarkably, the charge transfer is also accompanied by 
an orthogonality catastrophy. The analogy to the Kondo 
problem was recognized early with an exact formu- 
lation due to Matveev Q. For a single dot the Kondo 
type slow rearrangement of the electron states leads to a 
substantial downward renormalization of Ec [ 5 1 , as well 
as a smoothing of the conductance peaks |q, 7\. Addi- 
tional processes, such as the effect of higher order terms 
inelastic cotunneling ||, and mapping to the out of 
equilibrium Anderson model []l0|] were also analyzed. 

New effects arise when two such systems are allowed 
to interact. We argue that a genuinely new and robust 
quantum phase transition takes place in the coupled dot 
system when their mutual capacitance is varied. It is 
driven by a change of the degeneracy of the ground state. 
Interdot tunneling will be included perturbatively. It 
turns out to be a relevant operator, manifesting itself 
in an inverse power law temperature dependence of the 
total conductance at criticality. 

In the related system of two isotropically coupled 
Kondo impurities, a quantum phase transition was pre- 
dicted also, as a function of the interaction. The results 
of numerical renormalization group studies were con- 
firmed by conformal field theoretical methods and 
rationalized by phase shift arguments |i~3| ]. Unlike our 
case, the same Fermi sea of electrons interacts with both 
impurites and Particle-Hole symmetry is required to pro- 
tect the fixed point. This may render it harder to realize 
experimentally 

Let us start by considering a structure of two quan- 
tum dots, each coupled to their own leads. The lead- 
dot barriers are assumed to be narrow enough such that 
the tunneling is correctly modelled as a point contact. 
Furthermore we assume the presence of a strong enough 



magnetic field to achieve a fully spin-polarized electron 
gas. Thus the number of "flavors" , i.e. of the additional 
quantum numbers of transverse momenta and spin of the 
electrons is restricted to one. The dot is assumed to be 
relatively large, thus supporting a degenerate electron gas 
with small level spacing. This level spacing serves as a 
low energy cutoff, below which our scaling arguments do 
not hold. The Hamiltonian of one lead-dot system can 
then be written: 
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where e k is the energy of the electrons, J is the tunnel- 
ing amplitude and the indices a and take the values 
1, when referring to the lead and 2 when describing elec- 
trons in the dot. In a pseudospin notation the tunneling 
term is proportional to cr+g + a~g J|] . 

Next we recall that for small enough dots the Coulomb 
repulsion introduces an interaction term between dot- 
electrons, the scale of which is Ec = e 2 /2C, the charging 
energy, where C is the capacitance of the dot. Experi- 
mentally it is also possible to tune the overall potential of 
the system by a gate voltage Vq. The electrostatic energy 
of the dot can then be expressed as Eq = (Q — Qg) 2 /2C, 
where (essentially) Qg = CVq and Q is the charge on the 
dot. Tuning Qq beyond e/2 makes it energetically favor- 
able to transfer electrons across the barrier, giving rise to 
the well-known set of parabolae as the "band-structure" 
of the system. Transport becomes possible when the en- 
ergies of states with different number of electrons are 
degenerate. Thus the conductance shows sharp peaks as 
a function of Qg, with maxima at Qg/c = n + 1/2. In 
the vicinity of these degeneracy points the energies of the 
states with n and n + 1 electrons are much closer to each 
other than to any other state. It is then reasonable to 
truncate the Hubert space to two states. A second pseu- 
dospin of S = 1/2 can be introduced to represent this 
constraint on the allowed states. With this notation H 
assumes the Kondo type form: 
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Here A is the gap between the n and n + 1 electron 
states on the dot. The introduction of the two types 
of pseudospin operators allows a complete mapping of a 
single dot to the Kondo problem in a magnetic field A, as 
first realized in its entirety by Matveev Q . Note that the 
Kondo term is not spin-rotationally invariant, it contains 
only the spin flip terms. 

We proceed by including the interaction between the 
two dots caused by their mutual capacitance C m @,|l5| . 
The generated dot-dot coupling is proportional to til^r, 
where we introduced the L, R notation for the left and 
right dot respectively. Here ?il.r denote the charge of 
the left or right dot. In pseudospin notation >§£ — 
7t,l,r — 1/2. The mutual capacitance in pseudospin no- 
tation gives rise to an antiferromagnetic Ising type cou- 
pling: Hfg = J,S£S&, where I z ~ E Cm . The total 
Hamiltonian then takes the form: H = + + H^, 
describing two anisotropic Kondo impurities, coupled by 
an antiferromagnetic Ising term. 

The physical content of the model can be analyzed 
by generalizing the arguments of the theory of the two- 
impurity Kondo model At I z — we have two decoupled 
Kondo models. At T=0 in the magnetic language two 
independent isotropic Kondo singlets are formed with a 
"binding energy" ~ Tk, as the anisotropy of the Kondo 
coupling is known to be irrelevant around this fixed point. 
In the charge language, the electrons form strongly hy- 
bridized states between the lead and the dot. This hy- 
bridization manifests itself by the strong coupling Kondo 
phase shift, S = ir/2. The key observation is that the 
ground state is a singlet. At finite but small I z we fol- 
low Nozieres, who showed that all operators around this 
fixed point break the singlet and thus are of dimension 
two. In particular the dot-dot interaction involves virtual 
hopping operators to fourth order in the lead-dot hy- 
bridization amplitude. The corresponding diagrams con- 
tain a large number of fermionic operators, and thus are 
irrelevant. Alternatively the large number of fermion op- 
erators strongly confine the relevant phase space, leading 
to a positive exponent for the temperature dependence. 
This consideration again yields a vanishing contribution 
at T=0. 

In the opposite limit, I z — oo, the dot spins are aligned 
antifcrromagnetically. The (t, i) and Q, T) states are de- 
generate and form a doublet, which is independent of the 
conduction electrons. In the charge language these states 
consist of one extra electron being either on the left or 
on the right dot: (1,0) and (0,1). The energy of forc- 
ing on or taking away a dot-electron is ~ I z , and is thus 
prohibited in this limit. Let us recall that the interac- 
tion term between the lead electrons and the pseudospin 
contains only spin raising and lowering terms. In the 
allowed Hilbert space the matrix elements of this cou- 



pling are zero. Therefore the phase shift of the conduc- 
tion electrons vanishes. The degeneracy of the ground 
state extends to large but finite I z couplings as well, 
since the flip from one state to another again requires 
a high power of fermion operators and thus is irrelevant. 
In other words, whichever dot the electron resides on, 
it hybridizes only with its corresponding lead. Since the 
dot- dot coupling does not allow for charge transfer, these 
fluctuations remain confined to the dot and lead on the 
same side. Therefore the energies of the doublet's two 
states renormalize symmetrically for finite I z , and thus 
the degeneracy is preserved. 

To sum it up, the ground state is a singlet for small 
values of I z , but changes its symmetry to a doublet at 
large values of I z . This change cannot be continuous: the 
two regions are necessarily separated by a phase transi- 
tion. While this transition seems to be related to the 
case of two Kondo impurities, it is certainly different and 
its nature is as yet unexplored. In the remaining of this 
paper we assume that the transition is of second order. 

Here we pause to make connection to previous work by 
reviewing the band structure. The parabolae now have 
two indices, representing the charge states of the two 
dots. For I z — the (0, 0) and (1, 1) curves are touching 
Eq = 0, the latter displaced along the gate charge (Qg) 
axis by e. The (0, 1) and (1,0) curves are centered at 
Qg = e/2, and are also shifted upward such that they 
go through the intersection of the (0, 0) and (1, 1) curves. 
Exactly this degeneracy of states with different number 
of charges allows for transport across the dots, and thus 
gives rise to the conductance peak. If one now intro- 
duces the mutual capacitance I z , the upper parabolae 
are customarily shifted down by an amount ~ I z . This 
creates two degeneracy points at Qg ~ e/2(l ± I z /Ec)- 
Thus the original degeneracy of e.g. the (0,0) and (1,0) 
states, which allowed for the Kondo effect and underlied 
much of the above considerations, seems to have been 
destroyed. One might expect that whatever is left from 
the above picture, will be observable at the shifted de- 
generacy points. 

In contrast we predict that this new quantum critical 
point is observable, located precisely at the original value 
of Qg = e/2. The reason for this is that for small I z the 
Kondo energy scale Tk is bigger than I z . Therefore it is 
incorrect to construct a band structure first and then try 
to include the Kondo physics. Instead one has to start 
by accounting for the formation of the Kondo singlet, a 
deeply non-perturbative effect. The subsequent inclusion 
of I z then means only a small perturbation, similar to a 
fluctuating magnetic field. According to the above rea- 
soning such a field has a vanishing polarizing effect on the 
Kondo singlet unless its (Zeeman) energy is comparable 
to T K . Thus for I z < T K the (0,1) and (1,0) parabo- 
lae should not be viewed as shifted from their I z = 
location, and their degeneracy is preserved. A strongly 
analogous situation occurs in single dots: as shown first 
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in Ref. [|| , the effect of the Kondo type many body pro- 
cesses is to strongly collapse the band structure, sustain- 
ing their degeneracy up to some finite I z . On the other 
hand, for I z > Xk it is a reasonable starting point to ac- 
count for I z first, and then treat the Kondo coupling, i.e. 
the tunneling term as a perturbation. The two regimes 
are separated by the quantum critical point atij ~ Ik- 

Next we determine the total conductance of the two- 
dot structure. The interdot tunneling Hamiltonian takes 
the form: 

H tun =I±Y, C kh2 C k'R2 S L S R + h - a ( 3 ) 
kk' 

The pseudospin index 2 appears explicitly, as we are con- 
sidering dot-dot tunneling. This term breaks time re- 
versal symmetry, and is a relevant perturbation at the 
quantum critical point. One then expects that the low 
temperature behaviour of the renormalizcd tunneling at 
criticality exhibits a singularity: I±(T) ~ T~ 7 . The con- 
ducting path across the structure involves a lead-dot, 
dot-dot and dot- lead transitions. In the I± <C J limit 
the bottleneck, and thus the determining factor in the 
total conductance G is the dot-dot tunneling: 

G(I z =r z ,T)~(I ± (T)) 2 ~T~ 2 ^ . (4) 

This is only a crossover behaviour. As T is further low- 
ered, I± grows large and flows to an attractive fixed 
point, controlling its asymptotic behavior. Given its 
analogous structure, it is plausible that the dimension 
of I± is the same as that of the particle-hole symmetry 
breaking operator. However the actual value of 7 still 
needs to be determined p6[ |. 

What happens away from criticality? For I z < I c z the 
Kondo singlets inhibit the transport. At T = the bind- 
ing is complete, thus G(T — 0) = 0. Concentrating once 
again on the bottleneck dot-dot tunneling we compute 
the scaling dimensions of the involved operators. The 
fermion operators carry dimension 1/2, the spin rais- 
ing operator has dimension 1. The current operator is 
constructed from the [N, H] commutator. Collecting the 
terms the current-current correlator decays with the sixth 
power of time. Substituting this into the Kubo formula fi- 
nally yields G(T) ~ T 4 . The lead-dot process occurs via 
the Kondo coupling which scaled to its unitarity limit, 
thus it does not give rise to additional powers ot T. 

In the regime I z > I z electrons have to break an Ising 
bond. Thus at zero temperature again G(0) = 0, and at 
finite T the temperature dependence takes an activated 
form, G(T) ~ exp(-W/T), where W ~ I x . To sum it up, 
the conductance as a function of I z at zero temperature is 
zero nearly everywhere, and exhibits a resolution limited 
peak at I z = 1°. At finite but low temperatures the peak 
persists in G(T) as a function of I z . On the two sides 
of the peak G assumes non-zero values. These wings 
are asymmetric, with T dependent values. The different 
regimes are shown in Fig.l. 




FIG. 1. The conductance G as a function of T and I z . 

Finally we examine the effect of changing the gate volt- 
age, which tunes the gate charge away from the special 
point Qq = e/2, considered so far. In the Kondo lan- 
guage this gives a finite value to the magnetic field A. 
For I z < I z we utilize the same observation as before: 
for the unperturbed state is the Kondo singlet. Thus for 
A = 0, i.e. at Qq = e/2 transport is still impossible 
at T = 0. For < A < T K the singlet is somewhat 
polarized, and weak transport is possible. This mani- 
fests itself in two small-amplitude "shadow-bands" in a 
V shape determined by |A| = I z . This is the location of 
the crossing of the typically constructed "shifted parabo- 
lae" . An important transport channel in this region is 
co-tunneling, which only virtually breaks the Kondo sin- 
glet. For I z — I z the pronounced conductance peak of 
the quantum critical point is present at A = 0. This 
peak continues out to finite A, forming a parabola- like 
ridge, which smoothly connects to the usual split con- 
ductance peaks at |A| = I z for I z > I z - In this region 
I z is the larger of the energy scales and constructing the 
band structure first is appropriate. 

Constructing the picture from the large I z side, the 
magnetic field A is trying to induce a spin-flip transi- 
tion in the antiferromagnetic singlet. It is competing 
with the singlet binding energy, so the spin-flip can only 
occur when the binding energy equals the Zeeman en- 
ergy: |A| ~ I z , forming the usual V locus for the split 
peaks. Approaching the quantum critical point however 
the binding energy collapses to zero, hence the V becomes 
rounded, and closes at I z , as shown in Fig. 2. 

To summarize, the key experimental predictions of our 
calculations are: 

(i) If the gate voltage is fixed so that Qq = e/2, then by 
tuning I z a pronounced new conductance peak has to be 
observed at some critical value I z — I z ~ 7k- 

(ii) Staying at this point Qq — e/2, I z — 1°, the conduc- 
tance G(T) should exhibit a power-law singularity in its 
temperature dependence. 

(Hi) The amplitude of the split conductance peaks at 
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FIG. 2. The conductance peak as a function of I z and the 
gate voltage at T = 0. 



calculations are needed on this point. 

In sum we studied the system of two coupled quantum 
dots. We established the existence of an intriguing new 
quantum critical point. Several experimentally accessi- 
ble predictions were reached: a new conductance peak 
at Qq — e/2, an inverse power law dependence of the 
conductivity at this same point and a marked collapse 
of the split conductance peaks, when the experimental 
parameters are in the suitable range. 

We enjoyed illuminating discussions with I. Affleck, A. 
Georges B. Jones, A. Ludwig, A. Millis, T. Pohjola, A. 
Sengupta and G. Zarand. This work has been supported 
by NSF-DMR-95-28535 and the SFB 195 of the DFG. 



Q / e/2 should exhibit a marked collapse as a function 
of I z when I z approaches l c z from above. 

We know of no experimental observation of these pre- 
dictions yet. This maybe due to the fact that the above 
theory applies only under the following conditions: 

(i) Typical experiments |l7j and the corresponding the- 
ory (l4| considered the case of fixed I z , and described 
the evolution of the peak structure with tuning of I±, 
the dot-dot tunneling. The present theory addresses the 
case of fixed and small I±, and tuning with I z instead. 

(ii) The number of tunneling channels should be small. 
The above theory strictly applies only for the case of a 

single channel. This requires a narrow, long constriction 
between the leads and the dot, similar to the case con- 
sidered in We expect important changes when the 
number of flavors of the electrons is increased. Switch- 
ing off the magnetic field increases the number of chan- 
nels to two. It can be shown |l6| that the Ising term 
is marginal around the "decoupled" fixed point, leading 
to a line of fixed points which terminates at some inter- 
mediate value. The fixed point structure of the related 
two Kondo impurities model changes analogously. The 
Kondo singlet and the antiferromagnetic singlet phases 
remain intact, but the quantum critical fixed point ex- 
pands into a very unusual fixed point area, a whole region 
of the parameter space consisting of fixed points |l9[ |. If 
such a structure emerges in our case, then a broadened 
conductance peak will form as a function of I z at T = 0, 
and the finite temperature conductance should exhibit 
singular temperature dependence with I z dependent ex- 
ponents. The case of even larger number of channels has 
been investigated for single scatterers in relation to the 
physics of two level systems . It has been shown that 
a two dimensional subspace of the flavor indices emerges 
to dominate exponentially over the others in the course 
of scaling. Therefore we expect the basic picture of two 
distinct phases and a well defined quantum phase tran- 
sition in between to carry over, but obviously further 
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